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Study  of  the  Effect  of  Particle  Density  and  Particle  Size  in  /II 
a  Gas-Solid  Flow  on  the  Drag  Coefficient 
Zhou  Yanhuang,  Sun  Xingchang,  Yin  Hebao,  Chen  Xingquan  and 

Li  Guang 

[Abstract]  This  paper  shows  the  experimental  results  on  the 

drag  coefficient  between  the  gas  and  solid  pliase  within  the  range  of 

2  k 

Reynolds  number  Re  =  10  —10  ,  void  e  =  0.45—0.98  and  equivalent 
particle  diameter  dp  =  3— 11mm.  These  results  are  also  briefly 
analyzed  and  discussed. 

A  standard  drag  curve  has  already  been  given  for  the 

relationship  between  the  drag  coefficient  of  a  spherical  particle 

r  1  21 

and  Reynolds  number  ’  .  (See  Figure  6)  S.L.  Soo  et  al  be- 

[21 

lieved  that  this  result  could  be  applied  to  a  dilute  suspension  . 

In  a  packed  bed  loaded  close  to  its  maximum,  people  usually  use 

[2  3] 

the  result  obtained  by  Ergun  ’  .  In  interior  ballistics  of  a 

gun,  during  the  initial  combustion  period,  the  e  of  the  particle 
bed  is  around  0.5.  The  equivalent  diameter  of  the  particle,  dp, 

[A  ] 

ranges  from  several  to  a  dozen  millimeters.  E.B.  Fisher  ,  A.W. 
Horst^^^,  K.K.  Kuo^^^,  F.  Robbins  and  P.S.  Gough^^^  had  studied 
this  type  of  problem.  Their  joint  conclusion  was  to  continue 
using  Ergun's  results.  There  has  been  very  little  study  in  the 
drag  coefficient  in  a  scattered  suspension  with  e  ranging  from 
0.45—0.98,  which  corresponds  to  the  middle  and  late  stage  of 
combustion  in  interior  ballistics.  To  date,  we  only  find  the 
equation  proposed  by  K.K.  Kuo  which  is  based  on  the  Ergun  formula 
and  Anderson  formulate]; 


fl.7S 

(  «  <«.) 

1  1 

t. 

-f-J  (*.<  «  <ei) 

(1) 

,0.3 

(«, <«<!.•) 

where  is  the  friction  coefficient  calculated  based  on  the 
interfaced  area.  Cj  and  are  the  porosity  of  critical 
fluidization  and  the  maximum  porosity,  respectively.  Usually,  we 
choose  ei=  0.98.  This  equation  has  a  very  narrow  application 
range;  0.003<Re<2000. 

Therefore,  neither  the  standard  drag  curve  (Figure  6)  nor 
equation  (1)  can  satisfactorily  simulate  the  transport  phenomenon 
of  the  particles  in  interior  ballistics. 

In  order  to  verify  and  modify  equation  (1),  to  expand  its 
applicable  Reynolds  number  range,  and  to  understand  the  deviation 
of  the  drag  coefficient  of  the  kind  of  particles  of  interest  to 
us  from  Ergun's  results,  we  did  the  following  experiment  in  1981. 

Photographs  of  the  experimental  apparatus  and  some  of  the 
particles  used  in  the  experiment  are  shown  in  Figures  1  and  2. 

In  order  to  measure  the  drag  coefficients  of  loose  particles  at 
various  porosities,  c,  we  prepared  special  specimens.  Fine  nylon 
lines  were  used  to  suspend  particles  according  to  specific 
distance  and  direction  in  the  test  section  (See  Figure  3).  To 
the  extent  possible,  we  tried  to  make  sure  that  the  particles  are 
uniformly  distributed  so  that  they  are  isotropical.  The  method 
used  to  test  a  packed  bed  is  essentially  the  same  as  that  employed 
by  Fisher^^^  and  Horst^^^.  The  orientation  of  each  particle  is 
determined  by  the  random  process  used  for  filling  (See  Figure  4). 


2 


During  the  experiment,  the  pressure  drops  across  the  bed  at 
varying  height  were  measured.  In  addition,  parameters  such  as 
the  apparent  velocity  were  also  obtained. 


^Manuscript  received  on  March  13,  1984  and  revised  manuscript 
received  on  June  25,  1984.  This  paper  was  presented  at  the 
Second  Asian  Fluid  Dynamics  Conference  in  October  1983. 


Figure  1.  Experimental  Apparatus 


Figure  2.  Some  Particles  Used  in  the  Experiment 
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Figure  5.  Friction  Coefficient  vs.  Re  e/1-e 

1.  (1)  and  (2)  show  the  range  of  experiment  values 
obtained  by  the  authors 

2.  (3)  is  the  curve  recommended  by  the  authors 

3.  (4)  and  (5)  are  the  range  of  Horst's  experiment 

4.  (6)  is  the  result  of  the  Ergun  equation 
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Figure  6.  Drag  Coefficient  vs  Reynolds  Number  Re 
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1.  (1)  is  stand  drag  curve  for  spheres 

2.  (2)  is  Newton’s  internal  friction  law 

3.  (3)  is  the  curve  recommended  by  the  authors 

4.  (4)  is  the  curve  calculated  by  S.L.  Soo 

5.  (5)  and  (6)  show  the  range  of  experimental  values 
obtained  by  the  authors 
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Figure  7.  vs  Re  e/1-c  at  Various  Porosities 
1.  Ergun  equation 

Hence,  the  friction  coefficient  can  be  determined  by  the 
following  equation 


C. 


bP  /1-e  ■fPi 

-JJ-/—?  Zp 


(2) 


where  Ap  is  the  pressure  drop,  H  is  the  bed  height,  u  is  the 

s 

apparent  velocity  of  the  gas,  and  is  the  density  of  the  gas. 
The  drag  coefficient  Cjj  calculated  based  on  the  windward  area  of 
the  particle  is  related  to  as  the  following; 


(3) 


Sc 


We  measured  various  particles  of  11  different  geometric 
shapes  at  6  e  values.  After  considering  corrections  for  the 
nylon  line  and  the  mesh,  we  arrived  at  the  following  conclusions 
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When  e  =  0.45,  based  on  our  measured  results  and  on  the 

[2  3  51 

values  reported  in  the  literature  ’  ’  ,  the  relation  between 

and  Re  c/l-e  is  just  as  the  one  shown  in  Figure  5.  After  certain 
transformation,  the  relation  between  Cp  and  Re  is  shown  in  Figure  /14 
6.  At  various  e,  the  relation  between  and  Re  e/l-e  is  shown 
in  Figure  7.  Based  on  these  curves  we  can  see  that  our  results 
are  707o  higher  than  those  from  the  single  sphere  standard  drag 
curve  when  the  porosity  e  =  0.94~0.977.  At  maximum  packing,  e= 
0.45~0.47,  the  value  of  we  obtained  is  about  1.10,  for  less 
than  1.75. 

Based  on  our  measured  results,  the  following  relation  is 
introduced : 

C,.  (  e  <e,) 

C*.  (  1  _e,  ■  (e.<e<0.97) 


c,. 


0.45 


(0.97<  t<  1  ) 


where  is  defined  by  the  following  formula: 

[0.3l(lgffe)‘-2.55lg/?e+6.4  (J?e<20000) 


10  (7?^>20000) 


(5) 


The  value  obtained  from  this  formula  is  far  less  than  that  from 
the  Ergun  equation.  We  believe  that  this  discrepancy  is  due  to 
particle  size  difference. 

Due  to  the  presence  of  particle  clusters,  the  turbulent 
kinetic  energy  is  excessively  dissipated.  The  turbulence  is 
damped.  With  particles  of  the  same  size,  the  smaller  the 
porosity  the  smaller  the  gap  between  particles  becomes.  A  large 
scale  turbulence  is  being  further  divided  and  the  larger  the  drag 
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p 


coefficient  becomes.  In  a  suspension  of  the  same  porosity,  the 
smaller  the  particle  cluster  size  is,  the  smaller  the  average 
distance  between  particles  is,  and  the  larger  the  drag  becomes. 
Thus,  the  difference  between  our  results  and  those  of  Ergun  can 
be  fundamentally  understood.  Furthermore,  in  the  extremely 
dilute  case  with  e=  0.98,  we  chose  =  0.45  instead  of  0.30  as 
K.K.  Kuo  did  in  equation  (2).  This  was  based  on  the 
consideration  that  the  particles  we  used  were  cylindrical.  In 
addition,  the  authors  believe  that  even  when  e>0.98  the 
interference  among  particles  in  the  downstream  flow  field  is 
unavoidable . 
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A  New  Method  of  Rocket  Nozzle  Design 
Qiu  Guangshen 

[Abstract]  A  new  method  of  solid  rocket  nozzle  design  is 
introduced  in  this  paper.  This  method,  based  on  the  theory  of 
non-symmetric  flow  in  a  Lavol  nozzle,  is  first  focused  on  the 
reduction  of  aerodynamic  misalignment  and  then  on  efficiency.  I 
does  not  contradict  other  traditional  methods.  Instead,  it 
supplements  them.  This  design  method  will  help  reduce  the 
scatter  of  impact  of  free  flight  rockets. 

Table  of  Symbols 

<i)  -  potential  function; 

*  - 

critical  sound  speed; 

w  -  one-dimensional  flow  velocity  coefficient; 
h^  -  coefficient  of  <t)  in  flow  equation  as  represented  by 
dimensional  flow  parameters,  when  the  subscripts  i  = 

1,  2,  3  represent  xx,  rr,  and  rx,  respectively; 
f  -  small  perturbation; 

-  partial  derivative  of  with  respect  to  R.  The 
subscript  may  be  x,  e  or  x,  r.  Two  subscripts 
represent  a  second  order  partial  derivative; 

Jq,  zeroth  and  first  order  Bessel  function; 

the  first  root  of  J^(Bjj)=0  which  is  chosen  to  be  1.841; 

-  local  density; 

T-  -  axial  force  per  unit  area  at  the  outlet  end; 


-  transverse  cross-section  radius; 

Ma  -  one-dimensional  Mach  number  on  the  cross-section; 


dr  ^ /dx 


t 

k 


^Ma 


w 

^0 


a  function  of  x  alone; 

slope  of  the  bus  line  along  ox  direction,  which  is  tg  a 

in  the  linear  portion; 

radius  of  the  nozzle  throat; 

radius  of  curvature  of  the  throat; 

lateral  force  at  the  throat; 

specific  heat  ratio  of  the  combustion  gas; 

the  Planck-Myer  expansion  angle  on  a  cross-section, 

which  is  a  function  of  Ma; 

mass  of  the  propellant; 

total  mass  of  the  rocket. 


Manuscript  received  on  March  26,  1983,  revision  received  on  May 
23,  1984. 

1 .  Introduction  /41 

Because  free  flight  rockets  have  poor  density  of  impact, 
therefore,  finding  an  easy  and  feasible  method  to  improve  density 
has  been  a  major  subject  facing  the  free  flight  rocket 
technology.  As  far  as  rocket  design  is  concerned,  we  should 
start  with  reducing  the  eccentricity  of  thrust.  The  primary 
component  of  thrust  eccentricity  is  aerodynamic  eccentricity.  In 
recent  years  some  foreign  scholars  conducted  theoretical  and 
experimental  studies  on  the  non-symraetric  flow  in  a  Lavol  nozzle  to 
provide  the  basis  for  new  nozzle  design.  Normally,  the  principal 
objective  of  rocket  design  is  to  improve  efficiency.  The  major 
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nozzle  parameters,  to  some  extent,  are  freely  chosen  by  the 
designer.  The  new  design  method  based  on  the  above  theory, 
however,  is  focused  on  the  reduction  of  eccentricity,  as  well  as 
on  improvement  of  efficiency.  It  results  in  a  nozzle  whose 
theoretical  aerodynamic  eccentricity  is  zero.  The  nozzle  will 
facilitate  the  improvement  in  the  density  of  impact. 

2.  Brief  Introduction  to  the  Theory 

The  method  introduced  by  A.G.  Walters  is  to  treat  the  non- 
symmetric  three-dimensional  flow  in  the  nozzle  based  on  one¬ 
dimensional  flow  with  a  small  perturbation.  In  this  work,  the 
following  basic  assumptions  are  made.  The  expansion  section  of 
the  nozzle  is  axial  symmetric.  The  expansion  semi-angle  a  is  not 
too  large.  The  principal  lines  are  shown  in  Figure  2.1.  The  gas 
is  an  inviscid  gas.  The  flow  is  steady,  vortexless  and 
isotropic.  The  perturbation  is  small  and  plane  symmetric.  Due 
to  the  fact  that  a  is  not  too  large,  tgo  and  sin  a  can  be  treated 
as  small  first  order  quantities.  The  calculation  is  accurate  to 
the  first  order.  The  study  is  concentrated  on  the  flow  in  the 
expansion  section.  We  use  a  cylindrical  coordinate  where  the  ox- 
axis  coincides  with  the  axis  of  symmetry  of  the  expansion 
section.  The  origin  is  located  in  the  center  of  the  throat. 
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Figure  2.1  Nozzle  Shape 


First,  the  equation  of  the  flow  expressed  by  the  potential 
function  in  the  cylindrical  coordinate  is  linearized.  To  this 
end,  the  following  transformation  is  made 


(2.1) 

(2.2) 


Considering  the  fact  that  the  flow  near  the  throat  is  transonic, 
the  expansion  section  of  the  nozzle  is  divided  into  two  parts. 

The  intersect  of  the  curve  and  the  straight  line  portion  of  the 
nozzle  is  the  boundary;  i.e.,  x=Xo  .  Corresponding  to  X,  it  is  X# , 
After  neglectii.g  small  second  and  higher  order  terms,  we  get  the 
following  linear  equation 


*"*-x  * 


-( 


4>«k  + 


4 


•  <x<x. 


1 


(2.3) 


k-*  .-w  VT? 


Multiply  the  above  equation  by  cos  9  and  then  integrate  with 
respect  to  9  from  0-180®. 


«r‘+-^/ 

M 


(2.4) 


After  omitting  second  and  higher  order  terms,  we  get  the 
following  equation 

0  (2.5) 


Its  boundary  condition  is 


(//,)*U.,-  0 


Its  initial  condition  is 

(B) 

By  separation  of  variables  equation  (2.5)  can  be  solved  to  get 

m  m  \ 

F(R)  is  substituted  into  the  equation  for  A  to  determine  the  A 

n  n 

series.  A.  is  the  largest  one  and  the  rest  decrease  rapidly. 


Hence 


Q<X<X, 


(2.6) 


)a^XX~X,)-I^(9^X,y^^^(X^-X,)  X,<X  ^2.7) 


The  relation  between  L  and  is  determined  in  the 
following . 

The  thrust  is  moved  to  the  end  surface  of  the  nozzle  exit 
along  the  line  of  action  and  then  decomposed  into  two  components 
T  and  L  which  are  parallel  and  perpendicular  to  the  axis.  The 
moment  of  T  around  the  center  mass  G  of  the  entire  rocket 


(assuming  G  is  on  the  axis)  is  TZ,  which  is  denoted  as  as 
shown  in  Figure  2.1.  In  view  of  the  fact  that  the  inner  wall  of 
the  rocket  engine  is  most  frequently  a  body  of  rotation,  when  the 
nozzle  and  the  combustion  chamber  are  coaxial,  L  and  can  be 
expressed  in  terms  of  the  rate  of  change  of  the  momentum  and  that 
of  the  moment  of  momentum  of  the  gas  at  the  end  surface  of  the 
exit.  Then,  the  following  formula  can  be  obtained. 


Figure  2.2  Top  View  of  the  End  of  the  Nozzle  Exit 

(2.8) 

(2.9) 

where  the  symbols  relating  to  coordinates  and  angles  are  shown  in 
Figure  2.2.  Based  on  equations  (2.8)  and  (2.9),  by  omitting 
second  and  higher  terms,  we  can  derive  the  following  equation 

(2.10) 


w _  ^  \  2  Oj  dr, 

^  — 1ST  ~finF=\  -JT  “d* 


By  taking  a  small  element  of  the  expansion  section,  we  can  get 


another  formula 


(2.11)  /43 


Equations  (2.10)  and  (2.11)  are  the  fundamental  equations  to 
calculate  the  lateral  force  L  and  moment  C„. 


In  order  to  reduce  the  number  of  variables,  we  are  using  the 


following  relative  quantities 


TT 


Equations  (2.10)  and  (2.11)  can  be  re-written  in  the  following 


form 


/  •  f  u  \  dr,  r, 

^  A,y,  dx 

dC?  L* _ 

</x  r,  “  r,  dx  dx 


(2.10' ) 


(2.11') 


where 


Because  it  is  necessary  to  divide  the  expansion  section  into 
two  in  order  to  determine  ,  therefore,  it  is  also  required  to 
divide  it  into  two  parts  at  x,  (i.e.,  X#  )  when  we  determine  L* 
and  C^.  In  practice,  the  following  re-arranged  equations  can  be 
used  in  the  calculation. 


Asia  a 

Tf 


v'x*-x‘ 


0  <i<*. 


’  I  t*a 
I  +  A  X‘-x' 


0  <x<i. 


V  v-:  ••  v:-. : .  :v\- 


Near  the  throat,  i.e., 


at  X^Xa  , 


(D) 


1  +w;i 


X(k  +  1  ) 

The  relation  between  the  one-dimensional  flow  Ma  and  the  velocity 
coefficient  w  is 


_ 2h>* 

*  +  1  -(*  -  1  )u>* 


Far  away  from  the  throat,  i.e.,  Xjcx 

1 

dx  “  f,  1 

dM  _  tAfar  jrt  I 

d*  i/a*— 1  is  r, 

From  equations  (2.10')  and  (2.11')  we  get 


dCt 

dx 
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0  <5<X, 
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r  A/a*-2 
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0  <x<5. 
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Based  on  the  above  equations,  we  can  find  the  curves  showing 


L  and  vs.  x,  as  shown  in  Figure  2.3.  L  and  oscillate  with 
increasing  x.  The  amplitude  of  L  is  decreasing  while  that  of 
is  increasing.  The  phase  difference  between  the  two  is 
approximately  2/3lI. 


Figure  2.3  Lateral  Force  and  Moment  with  Varying  Expansion 
Section  Length 

We  verified  this  theory  experimentally  and  proved  that  it  is 
accurate.  For  instance,  a  comparison  is  made  with  respect  to  the 
x/r^  at  L=0.  The  result  is  shown  in  Table  2.1. 

Table  2.1  Comparison  of  Experimental  and  Theoretical  Values 


’  i 

1. 

«  ■  « 

1  •  «'  * 

•  «  • 

X 

1«.1( 

lO.Xl  ' 

•  »H 

1 

M.tr 

17. tl  1 

1.  experimental  value 

2.  theoretical  value 
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The  paper  by  A.G.  Walters  introduced  the  experimental  results  and 
arrived  at  the  same  conclusions. 

From  Figure  2.1  we  can  see  that,  if  we  assume  the  distance 
from  the  end  of  the  nozzle  outlet  to  the  center  of  mass  G  is  1 , 
then  the  moment  of  thrust  around  G  -  the  thrust  eccentricity 
moment  can  be  calculated  based  on  the  following  formula 

Af.=Z.;4-7.'=Z./+CT  (2.12)  /45 

From  the  correlations  in  the  L-x  and  C^-x  curves  we  can  see  that 
there  is  a  suitable  x,  which  can  be  chosen  as  the  length  of  the 
expansion  section,  to  realize  M^=0.  The  condition  is 


I  _  1  ^  __1_  ^2.13) 

-cr"“  /■  cj  *  i 

In  reality,  L  and  are  not  only  related  to  x,  but  also  to  other 
parameters  of  the  expansion  section  such  as  a.  and  R^,  as  well  as 
to  the  value  of  k  of  the  combustion  gas.  Therefore,  we  can 
choose  a  set  of  appropriate  parameters  to  meet  the  requirement 
=0 .  For  convenience,  we  suggest  that  this  type  of  nozzle  whose 
theoretical  thrust  eccentricity  moment  is  zero  should  be  called 
the  zero  eccentricity  moment  nozzle.  The  results  of  the  flight 
experiment  were  also  introduced  in  A.G.  Walters'  paper.  The 
presence  of  this  type  of  nozzle  was  proven. 

3.  Nozzle  Design 

The  principal  lines  shown  in  Figure  frequently  used  in 

free  flight  rockets.  The  advantage  is  ease  of  fabrication. 
However,  irregular  aerodynamic  misalignment  may  result,  which 
will  hurt  the  density  of  impact.  Because  of  the  sharp  edge  at 
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the  front  and  rear  end  of  the  throat  and  at  the  intersect  of  the 
contracting  cone  and  expansion  cone,  it  is  easy  to  be  eroded  by 
the  combustion  gas  to  form  irregular  arcs,  leading  to  irregular 
changes  of  L  and  C.^,.  Hence,  there  is  a  need  to  use  a  circular  arc 
at  the  interface  to  eliminate  erosion.  The  principal  lines  of 
the  nozzle  used  in  the  theoretical  derivation  are  shown  in  Figure 
3.1b  or  Figure  2.1.  There  is  no  cylindrical  section  at  the 
throat  which  obviously  does  not  favor  mass  production.  It  is 
more  appropriate  to  change  to  the  nozzle  line  (mainly  at  the 
throat  and  its  neighborhood)  shown  in  Figure  3.1c.  The  use  of 
the  R|  arc  as  a  transition  step  will  favor  the  reduction  of  the 
flow  perturbation  at  the  throat  and  the  improvement  of  the 
homogeneity  of  the  flow  at  the  throat  cross-section.  The 
expansion  section  of  this  type  of  nozzle  has  five  parameters:  r^, 
r^,  R^,  a  and  x.  As  far  as  the  principal  line  is  concerned,  r^ 
is  a  constant.  Only  three  of  the  remaining  four  are  independent. 
They  satisfy  the  following  relation: 

I +  -«■«)+ 1  (A) 

In  order  to  link  to  the  conventional  nozzle  design  methods,  we 
choose  5,  X  and  a  as  independent  variables. 
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Figure  3.1  Shape  Near  the  Throat 

In  the  initial  stage  of  rocket  design,  only  a  rough  shape  of 
the  nozzle  must  be  determined.  In  addition,  if  the  condition 
=0  is  used,  the  I  in  equation  (2.13)  is  still  an  unknown  or  is 
not  precisely  known.  Thus,  the  design  work  cannot  proceed.  In 
view  of  the  fact  that  for  most  free  flight  rockets  i  is  a  large 
number  and  the  calculation  shows  that  the  value  of  C,j,  is  very 
small,  the  requirement  Mj,=0  will  demand  that  L.l  and  arc 
equal  but  opposite  in  sign.  Thus,  we  can  determine  that  the 
point  which  satisfies  this  condition  must  near  L=0.  Hence, 
nozzle  design  can  be  carried  out  in  two  steps;  preliminary  design 
and  accuracy  connection. 

3.1  Preliminary  Design 

In  equation  (2.10),  there  is  a  factor  dr^/dx,  which  is  tgo, 
in  the  second  term  on  the  right.  It  is  a  small  quantity.  After 
neglecting  this  term,  based  on  the  fact  that  L  =  0,  we  get 


Therefore,  we  /46 


Obviously,  we  will  not  design  a  nozzle  with  x<Xo  . 
must  use  the  second  type  of  solution  for  1(3,  X).  Because  Aj 
( )  4  0 ,  then 


The  solution  obtained  is 


/.O.A-.)  J 


(3.1) 


In  the  formula  to  calculate  X,  there  are  h^ ,  h2  and  h^ .  After 
neglecting  second  and  higher  order  small  terms,  we  can  derive  the 
following 


X~X,+  — rr-t '' 
& 


(3.2) 


After  combining  (3.1)  and  (3.2)  we  get 


v(.Va)-/|S  tg-/f^(-Va*-T)-t6- 


(3.3) 

(3.4) 


The  initial  values  in  equations  (3.1)  through  (3.4)  can  be 
calculated  based  on  the  following  formulas: 


XMna 

(3.5) 

Ma,~  L 

rf*  +  i)-(*-i)ws 


Based  on  equations  (3.2)  to  (3.5),  we  can  get  a  family  of  C-a 
curves  using  X  as  the  parameter  at  fixed  k  and  n,  as  shown  in 
Figure  3.2.  Any  point  on  the  curve  represents  a  nozzle 
parameter.  This  parameter  satisfies  the  condition 

(Hi  )^_^  =0,  i.e.  L~0 

Based  on  Figure  3.2,  we  can  see  that  there  are  many  nozzle 
parameters  satisfying  the  above  condition.  Selecting  a  set  of 
values  from  these  combinations  of  parameters  must  also  be  based  on 
other  conditions,  such  as  improved  efficiency  and  structural 
requirements.  The  design  of  a  zero  thrust  eccentricity  nozzle  is 
to  introduce  a  constraint  of  nozzle  parameters  to  reduce  the 
aerodynamic  eccentricity  moment.  It  does  not  confuse  the 
conventional  design  methods. 

Due  to  the  fact  that  (H^  )___  =0  can  only  make  the  first  term 

r— ri 

on  the  right  of  the  L  formula  zero,  therefore,  the  lateral  force 
of  the  nozzle  from  Figure  3.2  should  be 

Mo‘-  1  r,  ix 

Using  the  value  of  x  corresponding  to  L=0  as  a  reference  point, 
because  <0 ,  therefore,  the  point  corresponding  to  must  be  on 
its  left.  Please  refer  to  Table  3.1  for  the  values  of  x  listed. 

If  M^=0  is  the  condition,  then  L  and  must  have  opposite 
signs.  Based  on  the  correlation  between  L-x  and  C^-x  curves  we 
can  see  that  this  point  must  also  be  located  to  the  left  of  the 
reference  point  (corresponding  to  L*0).  Hence,  the  nozzle 
prototype  determined  by  Figure  3.2  will  not  be  far  off  from  the  747 
zero  thrust  eccentricity  nozzle. 


Both  equations  (2.10)  and  (2.11)  are  obtained  by  the  same 


perturbation  factor.  In  reality,  there  are  numerous  factors 
capable  of  producing  perturbance.  Furthermore,  they  are  complex. 
Considering  the  small  perturbation  assumption,  the  effect  ol 
various  factors  is  independent.  The  combined  effect  is  the  sum. 
of  all  effects.  Therefore,  as  long  as  the  perturbation  satisfies 
the  initial  condition,  because  M  .=0,  then  EM  .=0.  Hence,  the 
theoretical  value  of  the  aerodynamic  eccentricity  moment  of  the 
nozzle  under  various  perturbations  is  still  zero. 

3.2  Accuracy  Correction 

The  prototype  nozzle  obtained  from  preliminary  design  does 
not  satisfy  the  condition  M^=0 .  As  the  design  work  proceeds, 
various  data  are  obtained  and  become  more  accurate.  After  the 


specific  conditions  are  established,  the  formulas  mentioned  above 
are  used  to  calculate  L*,  C^,  and  C^/L*.  The  set  of  parameters 
chosen  based  on  the  condition  C:|/L*=1/^  with  equation  (2.13)  are 
the  zero  eccentric  thrust  moment  parameters.  Prototype 


parameters  do  not  agree  with  this  fact.  We  must  rationally 
modify  it  to  satisfy  M^=0. 


Table  3.1  Comparison  of  Zero  Points 
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For  convenience,  we  can  refer  to  a  family  of  curves  or  a 
data  table.  Figure  3.3  shows  a  family  of  C-L  curves.  They  are 
calculated  from  the  C^/L*-x  curve.  The  variables  are  a,  X  and 
k.  If  C  of  the  nozzle  remains  unchanged,  is  known  (expressed 
in  terms  of  ^j^)»  and  both  intersect  at  A,  other  parameters  at  A 
can  be  determined  from  the  curve.  The  prototype  nozzle  can  then 
be  adjusted  based  on  these  parameters. 

Due  to  the  fact  that  the  theoretical  and  calculated  results 
are  very  close,  especially  the  difference  between  the  nozzle  line 
and  its  actual  shape  in  deriving  the  formula,  the  corrected 
nozzle  must  st?!!  be  experimentally  verified.  Because  the 
lateral  force  and  moment  are  not  too  large,  their  influencing 
factors  are  incidental.  Therefore,  the  magnitude  and  direction 
of  the  lateral  force  are  random.  This  special  feature  must  be 
considered  in  arranging  the  experiment. 
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3.3  Determination  of 

The  method  mentioned  above  is  essentially  using  the  baseline 
of  the  nozzle  to  control  the  line  of  thrust  action.  Once  the 
nozzle  parameters  are  determined,  the  intersect  between  the  line 
of  thrust  action  and  the  axis  of  symmetry  is  fired.  If  the 
intersect  coincides  with  the  center  of  mass  of  the  rocket  G,  then 
M^=0.  During  the  active  flight  period,  however,  the  charge  of 
the  rocket  is  continuously  consumed  and  exhausted.  The  center  of 
mass  G  is  not  fixed.  It  usually  moves  toward  the  head. 

Therefore,  it  is  not  possible  to  keep  M^=0.  Hence,  it  is 
necessary  to  determine  an  appropriate  I  value  to  satisfy  the 
design  requirements  for  zero  eccentric  thrust  moment  nozzles. 

Figure  3.4  shows  the  relevant  symbols.  G^  is  the  center  of 
mass  of  the  propellant.  is  considered  not  to  vary  with  time. 
Go  and  Gj^  are  the  center  of  mass  of  the  rocket  before  ignition 
and  after  the  propellant  is  combusted,  respectively.  Gj^  is 
chosen  to  be  located  between  G,  and  Gj^.  In  addition,  it  meets 
the  condition  that  the  misalignment  angle  of  the  total  thrust  is 


zero . 


Figure  3.4  Centers  of  Mass  and  Relevant  Symbols 
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Let  us  assume  that  the  instantaneous  center  of  mass  G(t)  is  at  a 
distance  a1  from  G,  and  Gj^  is  at  a  distance  from  G#  .  Based 
on  formulas  to  calculate  and  the  following  equation  can 

be  obtained  by  using  the  above  condition 


Qv^ 

rt  r: 


H  )df. 


(3.6) 


where  At  is  calculated  according  to  the  following  approximation 
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(3.7) 


The  upper  limit  of  integration  in  equation  (3.6)  can  be  selected 
in  two  ways.  One  is  to  choose  0.7'~0.8  of  the  time  corresponding 
to  the  critical  arc  length.  Or,  we  may  choose  t.  In  order  to 
use  the  value  of  exterior  ballistics,  the  former  upper  time 


limit  is  selected. 


can  also  be  calculated  by  the  following 


approximation 
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where  tj^  is  the  time  corresponding  to  Uj^ 
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The  curve  in  Figure  3.5  is  obtained  based  on  equation 
can  see  that  Ip  varies  very  slowly  beyond  u=4.5. 
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Figure  3.5  Ip  Vs.  Time 
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